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The recent demonstration of cooling of a macroscopic silicon nitride membrane based on dissipative 
coupling makes dissipatively coupled optomechanical systems promising candidates for squeezing. 

We theoretically show that such a system in a cavity on resonance can yield good squeezing which 
is comparable to that produced by dispersive coupling. We also report the squeezing resulting from 
the combined effects of dispersive and dissipative couplings; thus the device can be operated in one 
regime or the other. We derive the maximal frequency and quadrature angles needed to observe 
squeezing for given optomechanical coupling strengths. We also discuss the effects of temperature 
on squeezing. 
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I. INTRODUCTION 

The field of cavity optomechanics continues to regis¬ 
ter significant progress and a comprehensive review has 
recently appeared [l|. The important developments in¬ 
clude cooling of the mirror to its ground state , mode 
splitting i, quantum state engineering of the me¬ 
chanical motion p-Q, and electromagnetically induced 
transparency (EIT) and its various applications [iM3. 
More recently, strictly quantum effects like squeezing of 
the mirror [Tsl - E^ and the cavity field MS as well as 
generation of entangled photon pairs are receiv¬ 

ing considerable attention. Ponderomotive squeezing of 
light [27l - [3^ using an on-resonance driving laser is one 
of the most promising ways to generate squeezed light 
in cavity optomechanics. Safavi-Naeini et al. fabri¬ 
cated a micromechanical cavity resonator from a silicon 
microchip and they observed the fluctuation spectrum at 
a level (4.5 ± 0.2)% below the shot-noise limit despite 
highly excited thermal state of the mechanical resonator 
(10^ phonons). Purdy et al. [2^ placed a low-mass par¬ 
tially reflective membrane made of silicon nitride in the 
middle of an optical cavity and pushed the squeezing 
limit to 32% (1.7dB) by cooling the membrane to about 
ImK. Additional ways of producing optical squeezing in 
optomechanical systems have also been proposed. One 
example is use of a double-cavity optomechanical system 
to generate two-mode squeezed light [H, [s^. Another 
example [s^ is generation of quadrature squeezed light 
using the dissipative nature of the mechanical resonator 
in a single cavity driven by two differently detuned lasers. 
In a closely related subject, Lehnert and co-workers re¬ 
ported the experimental realization of entanglement be¬ 
tween cavity output photon-photon pairs [37| and en¬ 
tanglement between mechanical motion and microwave 
fields 

It should be noted that much of the work on cavity 
optomechanics uses dispersive coupling. However, there 
are a few studies for dissipative coupling p0l - l47l| — the 
intrinsic cavity lifetime depends on the mechanical mo¬ 
tion. A theoretical analysis of dissipative couplin g in 
cavity optomechanics was reported by Elste et al. j^ . 


They pointed out that the system gives rise to a remark¬ 
able quantum noise interference effect which leads to the 
Fano line shape in the backaction force noise spectra. Ex¬ 
perimentally, Li et al. [4l| for the first time reported dis¬ 
sipative coupling in a cavity optomechanics system that 
comprises a microdisk and a vibrating nanomechanical 
beam waveguide. Based on such a setup, Huang and 
Agarwal proposed a scheme to beat the standard 
quantum limit (SQL) by irradiation of squeezed light 
into the cavity. Hammerer and co-workers [i^, con¬ 
centrated on dissipative coupling by placing an optome¬ 
chanical membrane inside a Michelson-Sagnac interfer¬ 
ometer. This scheme is advantageous in the sense that 
the dissipative coupling is not due to internal dissipation, 
but the output photons are detectable. Weiss et al. 
presented a comprehensive study of dissipative coupling 
in both the weak- and strong-coupling limits, and they 
found the parameter regions for amplification of cool¬ 
ing as well as EIT and normal-mode splitting. Wu et 
al. [i^ experimentally reported the application of torque 
sensing by using dissipative optomechanical coupling in 
a photonic crystal split-beam nanocavity. Very recently, 
Sawadsky et al. (d^ demonstrated cooling starting from 
room temperature to I26mK based on the combined ef¬ 
fects of dissipative and dispersive coupling. This is quite 
a remarkable development where the couplings can be 
changed adding flexibility to the operation. Encouraged 
by the significant cooling in this experiment, we examine 
the optical squeezing that can be produced in a dissipa¬ 
tive optomechanical interaction. 

In this paper, we develop analytically the theory of 
ponderomotive squeezing in cavity optomechanics with 
dissipative coupling. We show that the squeezing mag¬ 
nitudes with dissipative coupling are comparable to 
those achieved using dispersive coupling. This squeez¬ 
ing scheme broadens the scope of the quantum study of 
nonlinear interaction in optomechanics. Our proposal is 
based on the parameters reported in [d^; however, it is 
not limited to this system and is applicable to any op¬ 
tomechanical systems that can provide combined inter¬ 
actions. This squeezing scheme works in the unresolved- 
sideband regime, which has advantages in its easier sys- 
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tern fabrication requirements. Moreover, this particular 
parameter regime makes it feasible for obtaining squeezed 
light with low frequency mechanical oscillators, although 
thermal phonons are still an issue. We show that the 
system can generate a 3dB squeezed field by use of rea¬ 
sonable driving laser powers when the thermal phonon 
occupancy is as large as 1.5 x 10® (the corresponding 
bath temperature T = IK). The effect of a higher bath 
temperature can be offset by increasing the driving laser 
power. As a by-product, our theory explains the new in¬ 
stability region for small pump laser red-detunings which 
was discovered in the experiment 

The structure of this paper is organized as follows: In 
Sec. El we introduce the Hamiltonian of the optome¬ 
chanical system with both dispersive and dissipative cou¬ 
plings, and find the input-output relation for the cavity 
field. In Sec. ED we provide the analysis of the squeezing 
effects under purely dissipative coupling. We compare it 
with the conventional dispersive squeezing and show that 
they both generate squeezed output with similar magni¬ 
tudes but in different quadratures. In Sec. IIYI we study 
the effects of the combined coupling on the squeezing and 
find the optimal quadrature angle for squeezing. We also 
study the effects of the mechanical mode at finite temper¬ 
ature. In Sec. El we analyze the effective detuning of the 
driving laser due to the change of cavity resonance fre¬ 
quency, and then show its effect on the squeezing spectra. 
We present our conclusion in Sec. ED 


II. MODEL 

We consider an optomechanical system in which a me¬ 
chanical oscillator (frequency ujm) is coupled to an elec¬ 
tromagnetic cavity. We model the cavity mode with 
the annihilation operator a and the mechanical oscilla¬ 
tor with the displacement x and momentum p, or with 
dimensionless operators Q = xjxzPF and P = (fixzPF)p 
where a:zPF = \/is the mechanical zero-point 
fluctuation. The mechanical displacements weakly mod¬ 
ulate the cavity resonance frequency uidQ) and damp¬ 
ing rate k{Q). We expand them to linear order to get 
uJc(Q) = ujc — QujQ and k{Q) = k — g^Q, where the dis¬ 
persive coupling constant gui = {dwc/dQ) and the dissi¬ 
pative coupling constant = (dn/dQ). In the general 
cases, the dispersive coupling is larger than the dissipa¬ 
tive coupling by a factor oi gi^/g^, = lOc/k^ 1. However, 
by placing a micro-membrane inside a Michelson-Sagnac 
interferometer, it has been shown that g^ and g^j can be 
made of the same order. 

When the optomechanical system is driven by a strong 
laser with frequency uji and power V, the Hamiltonian 
can be written, in the rotating frame, as 

H — Qi -f -t- P ) ligojix aQ 

-I-ill\/2 k(^[ a'^(f;-I-Oin) - i?.c.], (1) 


where £i 



and Oin represents the input vacuum 


noise. To proceed, we linearize the Hamiltonian following 
the standard procedure by writing a = Qg + ai, P = 
Pg + Pi and Q = Qg + Qi- The mean values of the 
steady state can be calculated as 


Qg 


y/^g£l 
Kg -I- iAg ’ 


Qs = ( — 




( 2 ) 


and Pg = 0. Under the effect of driving laser, the me¬ 
chanical oscillator displacement Qg modulates both the 
cavity resonance frequency and the decay rate. Hence we 
define Ag = (wc — g^Qs) — w/ as the driving laser detun¬ 
ing from the effective cavity resonance frequency; and we 
define Kg = k — g^Qg as the effective cavity decay rate. 
Both Ag and Kg depend on the power of the driving laser. 
However, by tuning the driving laser frequency w/, one 
can always cause it to be on resonance with the effective 
cavity frequency, i.e. Ag = 0. Under this condition, the 
effective cavity decay rate is determined by the quadratic 
equation — KKg + 2£fgi^gi^/L>Jm, = 0. In the typical op¬ 
tomechanical systems, the term 2£fg^gi^/ijJm is negligible 
compared to k and hence Kg = k. For example with the 
parameters reported in [i^, 2£fgi^g^/uJrn < k/ 10® when 
the driving power is below lOmW. 

Then the linearized Hamiltonian takes the form H = 
Hq -f Pint T Pdamp and 


Ho = hAga\ai + -huJm{Q\ + Pi^), 

Pint = - /iG«%^Qi, (3) 




2y/2i 


Pdamp = ifiA/2^(a|ain - a|'„ai) - 


where = V2asga;,K is the driving field enhanced 

dispersive (dissipative) coupling constant. The form of 
the Hamiltonian (|31) suggests that it is more intuitive to 
write the cavity field in terms of its quadratures: X = 
(oi -I- a\)l\/2, Y = (ai — a\)/{y/2i), and [Ar,y] = i. 
In this paper, we are interested in generating squeezing 
light in the output field. Under the effect of dissipative 
coupling, the standard input-output relation 

K - gK.Qs) 

= y/^ai - (4) 

v2Kg 

since g^ <C Kg. In terms of quadratures, the input- 
output relations are Ajn + Aout ~ \/2j^X — Qi and 

Yin + Yout ~ \/2KgY. Hereafter, we first focus on the 
on resonance driving scenario (Ag = 0) and then discuss 
the squeezing effect with detuned driving by relaxing this 
condition. When Ag = 0, the coupling strength is 
real. The dynamics of the system can be described using 


T ^out ~ \/2(i^ g^Qs)^i 
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FIG. 1. The input quantum noises and their coupling rela¬ 
tions among different quadratures {X, Y) of the cavity field 
and mechanical mode (Qi). The dashed arrows show the 
noise input and coupling due to dissipative coupling G^. 

the quantum Langevin equations 

—Qi + —Qi + u:mQi = G^x + + c, 

OJm ^ V 

( 5 ) 

X = -KsX + (6) 

Y = -KsY + G^Qi + (7) 

Here, ^ models the Brownian noise acting on the me¬ 


chanical oscillator, and it obeys = 7m(2hth + 

l)(5(t — t'), where hth is the mean phonon occupation 
number. The correlations for the vacuum field are 
2Ks{Xin{t)Xin{t')) = 2Ks{Yin{t)Yin{t')) = K.s5{t-t'). In 
the unresolved-sideband limit Ks ^ <-^m ^ 7m, hence the 
vacuum noise dominates over the Brownian mechanical 
noise at low hth- 

We illustrate the coupling relations of the quantum 
noises in the optomechanical system, in Fig. [TJ The field 
quadratures are subjected to the vacuum input noise X-^ 
and Ijn- More importantly, we notice that, due to the 
dissipative coupling G^, the input vacuum noise Ijn is 
also coupled directly to the mechanical motion Qi. At 
the same time, the form of the interaction Hamiltonian 
shows that Qi interacts with the different cavity quadra¬ 
tures at the rates G^, and G^. Therefore, Ijn is fed into 
the system through two paths: (i) it directly couples to 
the cavity field; and (ii) it couples to the mechanical mo¬ 
tion Qi dissipatively and then the optomechanical inter¬ 
action transfers the noise to the cavity field. These two 
paths interfere in a coherent manner and lead to the Fano 
resonance in the cavity field spectrum. 

We calculate the output field by combining Eqs. 0-© 
after taking the Fourier transform, and find 


/1 G^G^ Ks+ioj ^ ^ 

I-- ~^m) — 

X Us — 

,1 GoiGk Ks+iuj . 

( : Ji^out : ^inj 

X Us — lUJ Ks — 


iioGfiGi^ 

{ks - iuj) 
2^sGl 
(ks - ito) 


where % = ojm/(<^m — — iujjm) is the mechanical sus¬ 

ceptibility. Equations® and ® describe how the input 
quantum noises add to the quantum fluctuation of the 
output fields. Without optomechanical interactions, the 
output field preserves the input field fluctuations, i.e., 
{XLt) = (F^ut)- 0116 increases the optomechanical in¬ 

teraction strengths G^i and G^, the noises are distributed 
in a nonlinear manner. The quantum squeezed states are 
generated when the variance is lower than the that of 
the coherent state, i.e., Sg = (Zg) < 1/2 for a specific 
quadrature Zgo^t = Xout cos 9 + F^ut sin 0. 


III. SQUEEZING WITH PURELY DISSIPATIVE 
COUPLING 

The phenomenon of ponderomotive squeezing with 
purely dissipative coupling can be obtained by setting 
the dispersive coupling strength G^j = 0 and = 0, so 

that Fout — Fjn and Xout = FjnT mechan¬ 

ical noise. The vacuum input Fj^ is coupled, not only 
to F^ut, but also to Xout via the mediated mechanical 
mode Qi scaled by the mechanical suspectibility x 


Win - 


2Ks{Ks — fw)' 


rFjn + 


iujGk 


V I i^G^Gu: 
7-d^iii“r/ . \o-^in“r 

(aCs - 


\/2ns{Ks — ioj) 
G^'\/2 ks 




Ks — liO 




( 8 ) 

( 9 ) 


dissipative coupling strength G^. When one measures 
the field ^gout = Xout cos 9 + Fj,ut sin 0 at 0 yf 0 ° or 90 °, 
Fout interferes partially with Xout since x(a;) is generally 
complex. The interference leads to squeezed quantum 
noise. The output squeezing spectrum is 

•S'diss - ^ + (^IxPrdiss cos^ 9 - Rex sin 26>), (10) 

where Fdiss = G^ux^j (dtcf) -|-7m(2hth +1) is the effective 
mechanical damping rate. By optimizing 9 and x(<^) ive 
obtain the optimal squeezing magnitude 


oopt _ 

*^diss ~ 


_7m(^lth + 1)_ 

Glxil,l{^K?^ + 27„^(fith + 1) ■ 


( 11 ) 


The squeezing magnitude can be enhanced by a large 
effective dissipative optomechanical coupling strength 
G^w^/ (iK^Xm) and a low mean phonon occupancy num¬ 
ber hth- The optimal squeeze d quadrature angle lies at 
^ ^/Glu}l,/{2Kl^rn), and 6»°P* approaches to 
90 ° with a large dissipative coupling strength G^ - From 
the above analysis, we can see that the ponderomotive 
squeezing relies solely on the interference of two paths of 
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Xin- One needs to suppress the input noises and ^ 
by choosing a quadrature angle close to 90°. The 
output field shows anti-squeezing at a; = ujm when 9 0. 

To illustrate the squeezing effect, we plot the output field 
spectra at different quadratures in Figs. [2ja) and [2jb) 
by numerically solving the quantum Langevin equations 
m-m- We use the parameters provided by the experi¬ 
ment reported in |48l |. and the specific values are given 
in the caption of Fig. [H At the angle the output 

spectrum [as shown in (b)] is characterized by a large 
squeezing of ~ 20dB at frequency uj ^ uJm — 27r x 2Hz 
and anti-squeezing at w = w^,. 



Frequency {u;^ -lj)/ (kHz) Frequency -cj)/ (kHz) 


and x(w) we obtain the optimal squeezing magnitude 


QOpt _ 

‘^disp ~ 


7m (nth + 1) 

+ 27m(fith + 1) 


(13) 


This result is identical to the one derived in [3(1 and has 
been experimentally demonstrated in . The opti¬ 
mal output frequency is (w - = FdispTm/S -f 7m/4, 

which increases with coupling strength G^. We plot the 
output spectra of dispersive squeezing in Figs. [IJc) and 
Hd), as a comparison with the dissipative squeezing in 
Figs. Ha) and Hb) ■ The optimal squeezing spectrum 
has a quadrature angle close to 0. The optimal squeez¬ 
ing magnitude is shown as ^ 40dB, which agrees with 
Eq. (fT^ . We observe similar output squeezed spectra, 
but the optimal squeezing magnitude is larger than in 

Figs. Ha) andHb). 

Physically both the dispersive coupling and the dissi¬ 
pative coupling generate optical squeezing in a similar 
manner, in the sense that they couple the input noise 
from one quadrature coherently to the other quadrature. 
Thus the input vacuum noise couples to the optomechan¬ 
ical system via two paths, as shown in Fig. [TJ These 
two paths interfere and lead to squeezing. The optimal 
squeezing exists at different quadrature angles due to the 
fact that Gui couples noise from A to T and G„ couples 
noise from T to A via the mechanical mode. 


IV. SQUEEZING WITH COMBINED EFFECTS 
OF DISSIPATIVE AND DISPERSIVE COUPLING 


FIG. 2. (Color online) Comparison of squeezing spectra with 
purely dissipative coupling (a) and (b) and purely dispersive 
coupling (c) and (d): The regions inside the black contours 
in the density plots (a) and (c) show a 3dB squeezing region 
and the blue horizontal lines show the optimal quadrature 
which are plotted in (b) and (d), respectively. The dissi¬ 
pative coupling strength is Gk = 2n x 150kHz with driv¬ 
ing laser power V 3.5W; the dispersive coupling strength 
is Gui = 2n X 75kHz with driving laser power V ~ 40mW. 
Other parameters are Ks = 2n x 1.5MHz, LJrn = 27r x 136kHz, 
7m = Stt X 0.23Hz, As = 0, and hth = 0. 


In the other limit when dispersive coupling solely gov¬ 
erns the optomechanical interaction, j.e., = 0, Eqs. 

dH) and (H reduce to Aout = Xi^ and Tout = -|- 

{xGl/Ks)X,n+ mechanical noise. This is the conven¬ 
tional ponderomotive squeezing scheme. It shares a sim¬ 
ilar noise transformation with that we discussed above. 
Hence we are able to observe a similar squeezing phe¬ 
nomenon, but the optimal squeezed quadrature is around 
tan6»[[Pp ^ \/Ks 7 m/( 2 G 2 ), and approaches 0 with 
a large dispersive coupling strength Gi^. The output 
squeezing spectrum is 

S'disp - + —(2|xprdispSin2 6»-bRexsin26»), (12) 


where Fdisp = G^/ks + 7m(2nth + !)• By optimizing 9 


In the previous section, we studied squeezing phenom¬ 
ena with purely dispersive coupling or dissipative cou¬ 
pling. One natural question is whether the combined 
effect of these two coupling regimes could enhance the 
squeezing. We next study the generation of squeezed 
state in the presence of both coupling regimes G^j and 
Gk- The output squeezing spectrum is 


*S^comb — 


1 


2|G|2 

-b 




X 

F' 


Re 


F,. = 


m- 

iGt . 4G2 


+ 




F,^ sin^ 0 cos^ 9 


2|£i|2 

sin 20, 


(14) 


■7?7i(2u.th 4” 1); 



4, ,4 


4k® 

2GS 


- 3 —7m(2?^th + 1), 

GW 

2k3 ’ 


where D = 1/x ~ GuiGk,/{ks — ioj). The optimal squeez¬ 
ing quadrature angle tan20[[y^j. = (F^j — r„)/(2r') and 
the squeezing magnitude can be found to a rough approx¬ 
imation 


oopt ,-s^ 

‘^comb 9 

^ I u 


r'V2 


+ r« + 2v'r,,r, + r' 27 ^’ 


( 15 ) 
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FIG. 3. (Color online) The density plots (a) and the optimal 
squeezing quadrature (b) of the output field spectra with com¬ 
bined dispersive and dissipative couplings. The regions inside 
the black contours in (a) show the 3dB squeezing region and 
the blue horizontal line shows the optimal quadrature which is 
plotted in (b). The coupling strengthes are Gui = 2n x 75kHz 
and Gk = 27r X 15kHz with driving laser power V ~ 40mW. 
Other parameters are identical to those used in Fig. 0 

In Fig. |3Ka), we plot the output spectra at different 
quadratures when the optomechanical system is subject 
to both dispersive and dissipative couplings. We set the 
coupling strengths such that = 5Gk in accordance 
with the experiment parameters in [48j | . The density plot 
resembles the main feature of ponderomotive squeezing 
with purely G^j or G ^, except for a trivial quadrature dif¬ 
ference. However, there are distinctions. The frequency 
bandwidth of the squeezing spectra increases at large 
quadrature angle and shrinks at lower quadrature angles. 
This is particularly advantageous in practice, since one 
usually focuses on a specific quadrature and hence one 
can make use of the larger bandwidth of the squeezed 
spectra. 

In the optomechanical ponderomotive squeezing pro¬ 
cess, the mechanical element functions as an active me¬ 
diating element and it provides coherent coupling be¬ 
tween two field quadratures. At the same time, it is 
subject to the environmental Brownian noise which is 
incoherent with the cavity field. In the reported pondero¬ 
motive squeezing experiments with purely dispersive op¬ 
tomechanical coupling, the environment temperature sets 
the limit of the squeezing magnitudes: Safavi-Naeini et 
al. re ported 0.2dB squeezing at hth 10"* and Purdy 
et al. |29l| pushed the squeezing magnitude to 1.7dB with 
a lower thermal phonon occupancy hth = 47. 

We now compare the effect of the thermal phonons 
on squeezing with different optomechanical couplings. 
Equations. (US, (m and (13 indicate that the out¬ 
put quadrature variance increases approximately pro¬ 
portionally to hth at large coupling rates. Comparing 
Eqs. (fTTl) and (fT^ . we find that optomechanical systems 
with purely dissipative coupling (Gk) or purely disper¬ 
sive coupling (Gcj) can generate squeezed field of similar 
squeezing magnitude if G^ujm = V^GojK-s. In Fig. |4l we 
illustrate the effects of the mean thermal phonon number 
on the optimal squeezing magnitude under different cou¬ 
pling regimes. The curves show that the squeezing mag¬ 


FIG. 4. (Color online) The effects of the mean thermal 
phonon occupation hth on the optimal squeezing magnitudes 
with different couplings. The optimal squeezing magnitudes 
are very similar for finite hth; hence the three curves overlap. 


nitudes decreases with large thermal phonon occupancy 
hth • Even when the thermal phonon number is as high as 
hth = 1000, the system yields about 15dB squeezing with 
combination optomechanical couplings at 7^ = 40mW. If 
we increase the driving laser power to V = 280mW, the 
squeezing magnitude increases to 30dB. Note that this 
phonon number is however difficult to achieve with low 
mechanical frequency Wm since hth is inversely propor¬ 
tional to LOm- For example, the system has to be pre¬ 
cooled down to T ~ 6.5mK in order to get hth = 1000. 
On the other hand at high bath temperature, a large 
squeezing magnitude requires an increase in the coupling 
strength, which can be achieved by increasing the pump 
power. If the bath temperature increases to T = IK, 
the corresponding thermal phonon number increases to 
hth 1-5 X 10®. One needs to increase the driving laser 
power toV^ 280mW in order to get 3dB squeezing. For 
this power the system is still in the stable region. 

Note however that the pump power cannot increase 
infinitely as too strong a pump laser leads to instability of 
the system dynamics. We discuss the stability condition 
in detail using the Routh-Hurwitz criterion (4^ in the 
Appendix. For example, our linearization method breaks 
down and the system settles into instability when the 
laser power reaches V ^ 290mW for the parameters given 
above and the laser frequency set as w; = wq — ujm- At 
this power, the coupling strengths are Gi^ = 27r x 200kHz 
and Gk = 27 r x 40kHz. A lower driving laser frequency 
allows for a higher critical pump power. Our analysis 
also reveals an unstable region when the effective driving 
laser detuning Ag is a small negative value. This explains 
the special instability region discovered in [4^ . 


V. SQUEEZING WITH A FIXED FREQUENCY 
DRIVING LASER 

Sawadsky et al. [i^ demonstrated a strong cooling ef¬ 
fect in an optomechanical system with both dissipative 
and dispersive coupling interactions. The experimental 
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results agree remarkably well with the theoretical calcu¬ 
lation. In the experiment, the authors fixed the driving 
laser frequency w; on resonance with the empty cavity 
resonance frequency Wc • When the driving laser power in- 



FIG. 5. (Color online) The change of the effective detuning 
and mean cavity photon number as the driving laser power 
increases from 0 to 200mW. Other parameters are identical 
to those used in Fig. [3] 


VI. CONCLUSION 

In conclusion, we investigated the generation of 
quadrature squeezed states with dissipative coupling op¬ 
tomechanical interactions. Our results show that the dis¬ 
sipative coupling interaction is able to generate strong 
squeezed vacuum states. The squeezing magnitude de¬ 
pends on the coupling strengths and the mean phonon 
Squeezing Spectra of Combined Coupiing 




FIG. 6. (Color online) The squeezing spectra in an optome¬ 
chanical system with both couplings combined with = 
5Gk = 27r X 75kHz and V 40mW. Different ranges of fre¬ 
quencies and quadrature angles are shown in (a) and (b). 
Other parameters are identical to those used in Figs. [2] and 
[3l The regions between the black contours have over 3dB 
squeezing. 


creases, the effective cavity resonance frequency changes 
due to the displacement of the mechanical membrane and 
this leads to an effective detuning of the driving laser. 
In this section, we analyze the squeezing phenomena in 
an optomechanical system driven by a laser with fixed 
frequency w; = Wc- Under this condition, the effective 
detuning and effective cavity decay rate Ks can be 
determined by solving the nonlinear equation set ([2]). We 
use the parameters reported in [i^. The solution to ([2]) 
shows that Kg k when the driving laser power Ei is be¬ 
low 250mW. However, the effective driving laser detuning 
As increases linearly from 0 to a value close to —uJm, as 
shown in Fig.[5l The cavity mean photon number |as| is 
also displayed in Fig. [5l When the driving laser power is 
set as 40mW, the effective detuning Ag = 27r x 20kHz. 
The corresponding coupling strengths remain the values 
Gcj = 27 r X 75kHz, which are similar to the ones used 
in Figs. [2] and [3l We show the squeezing spectra with 
combined coupling interaction in Fig. [6] at zero tempera¬ 
ture. The optimal squeezing magnitude reaches close to 
40dB. We find large regions with over 3dB squeezing in 
the spectrum, as illustrated between the thick black 3dB 
contour lines. We observe large regions of squeezing over 
lOdB and even squeezing over 20dB in Fig.[6l The result 
is very similar to the one in Fig. [31(a), and even the effect 
of temperature is similar so it is not discussed here. 


occupancy due to the mechanical noise. When the dissi¬ 
pative and dispersive coupling strengths are similar, they 
both generate comparable squeezing magnitude. This 
scheme works in the unresolved-sideband limit which en¬ 
ables its application in low-frequency mechanical oscil¬ 
lators. In potential experimental realizations, one chal¬ 
lenge would be the large thermal noise introduced by 
the large phonon number with low mechanical frequency. 
Large squeezing magnitudes require one to precool the 
system using a dilution refrigerator. The large thermal 
noise can also be offset by increasing the pump power by 
the same order of magnitude as for nth- 


APPENDIX: STABILITY CRITERION 


While we follow the standard linearization procedure 
in solving the nonlinear Hamiltonian, we must make sure 
the stability of the system dynamics for our chosen pa¬ 
rameters. We investigate the dynamics of the system 
using the quantum Langevin equation 


d4'(t)/dt = ri4'(t) -b (16) 


with 'l'(t) = (X,Y,Qi,Pi)^ for the system operators. 


'l'in(i) — (V^^sAlin, 0, ^-b 


ItxiG. 

\/2 ks 


Xs„- 


ReG;, 

\/2ks 


Y r 

in^ 


for the input noises, and 
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( Ks 

0 0 

y Re[Gcj] Im[Ga;] + 


Re[G,] - Im[G^] - :^Si 
Iiii[Gk] + Re[Gtj] 

0 


0 

-7m / 


(17) 


Note Guj,k = system is stable if all the eigen¬ 

values of the matrix [M have negative real parts. This can 


be examined by applying the Routh-Hurwitz criterion to 
the polynomial of its eigenvalues. The Routh-Hurwitz 
coefficients are 


hi = 2 Ks + 7m, 
/l2 = 2jrnK‘s + 


hs — 7m (A^ + Kg) -|- 2KsU)^ + sf 


V^^kAs 4^q;Ks 


A2 -I- kI 


h4 = <(A2 + K^)-Agffg2^ + 


2\/2£fuJn 


2 ks kI + A^ 


+ Asg^){gf,As 


V2g^Ks). 


(18) 


Routh-Hurwitz criterion reads that the system is stable 
if and only if all the coefficients hi > 0 for i = 1,2,3,4 
and the determinants of all of the Hurwitz matrices are 
positive, i.e. 


Note that when —>■ 0, these conditions reduce to the 
stability condition for the optomechanical system with 
purely dispersive coupling g^^. 


hs > 0 , 


h4 > 0 , 
hih2 — ha > 0, 
^1/12^3 — h^hi — hg > 0. 
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